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ON LAGRANGE THEORY OF SHELLS OF REVOLUTION 


VLADIMIR V. ELISEEV, ALEXEY S. SKOBLIKOV 

Abstract. The new linear theory of elastic shells is presented in this paper. This theory 
is free from various logical imperfections, that may be found in the approaches of earlier 
researchers. On the base of this theory the equations of shells of revolution are built. The 
equations of cylindrical shell are compared by a number of properties with the equations of 
other researchers and with elasticity equations of 3-dimensional tube. 

Important quality decisions are stated. 


1. Introduction 

Tensor analysis in the theory of elasticity was widely used by A. I. Lurye (for example, [8]). 
The direct tensor analysis prevails: tensors are considered as invariant objects, component 
form is used rarely, only when necessary. This approach makes analytical calculations more 
comprehensible, letting us pay more attention to the sense of formulas, not to algebra. 

Recently great success was made in the theory of elastic shells. First approaches were 
made by Berdichevsky [2], then direct tensor algebra was involved by V. V. Eliseev [5, 234], 
Now it seems to be logically and mathematically perfect. It was not published in English 
yet, therefore we will first present this new approach, based on direct tensor analysis and 
principle of possible dislocations (Lagrange principle). 

Having general theory of elastic shells we can build simplified forms of the theory. Here we 
derive the equations of shells of revolution. The equations of cylindrical shell are compared 
by a number of properties with the equations, derived by other researchers ([3], [4], [14], [7], 
[6]) and with 3-dimensional model of the tube. 

Important quality decisions are stated. 

2. General Theory 

The bases of the shell theory are considered to be well known [7], [9], but its forming is not 
finished yet. Traditional approaches yield to asymptotic decomposition of 3D problem [7], 
[2]; nevertheless the presentation of the shell as material surface doesn’t lose its significance 
[1], [11], [15]. In this study we consider linearly elastic shells as surfaces of particles with 
five degrees of freedom — material normals. Rooting in methods of analytical mechanics 
and representation of stresses as Lagrange multipliers [12], it’s easy to derive all equations of 
Kirchgoff-Love type theory. The results slightly differ from the already known. 

2.1. Geometry of material surface. Radius-vector is the function of two coordinates 
r[q a ). Let’s introduce vector basis in the tangent plane r a = = d a r, normal ort is 

n = . reciprocal basis is r a ( r a ■ rp = <5g), Hamilton operator V = r a d a , and metric 

tensors are a = Vr = r a r a and b = —Vn = — r a n a , n a = d a n. Derivative formulas 

^a —0 = T.oi 0 ^a 0 \T. T b a pn, n a b a pj_ 

T q / 3 a = - (d a a/ 3 \ + dga a \ — d x a a p ) 
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show that the form of the surface is totally defined by the dependency of covariant components 
a a p and b a p from coordinates [7]. 

It’s important to define precisely the degrees of freedom of the particles of material surface. 
The particles of momentless shell have only three degrees — translation. In the Cosserat 
model particles are solid bodies with six degrees. But classical views of Kirhgoff-Love imply 
another model: particles are material normals with five degrees. Vector of small rotation 
6 in this case lies in the tangent plane, like all moments. Small change of normal during 
deformation and work are 


n = 6 x n = 0, m ■ 9 = m' • (j), m' = m x n 

(it’s more convenient to work with 0 rather than with 9). 

Vector of displacement u — is the increment of radius-vector (f). Its connection with 
rotation follows from orthogonality 

(1) r a ■ n — 0 =>■ u Q ■ n + r a ■ 0 = 0 =>■ 0 = —Vw • n. 

The form of the considered material surface is defined by covariant components a a p and 
b a 0 , its small deformations are characterized by tensors 

£ = ]^a a 0 r a r p , % = b a 0 r a r' 3 . 

Evaluating the components’ increments we get 

daft = Ha ■ Hi3 ^ bj a0 = u a ■ Vp + Ha ' Htf} ^ 

=>- e — - (Vu ■ a + a ■ Vu T ) = (Vm s )ii 

( 2 ) = 2 - = = - 

b Q0 — n a ■ t 0 b a p — 0 a • rp n a ■ up 

=>• x = — (V0) + b ■ Vu T 

(sign (.. ,)|| means belonging to the tangent plane or can also be expressed as v ■ a for vectors 
and a ■ T ■ a for tensors of rank 2). 

Strain tensor e, defining the change of lengths and angles on the surface, takes place in all 
presentations of the shell theory, but tensor x in the form (2) is not so well-known. Let’s 
note that both tensors are symmetric and belong to the tangent plane. Due to the relation 
(1) tensor of bending and rotation x can be presented like: 

x = (VVu • n) || 

(it resembles classic theory of plates). 

2.2. Virtual work principle and its consequences. It’s possible to derive all system 
of equations of elastic system from this differential variational principle (except equations 
caused by material structure). For the considered model of the shell the following form of 
the principle is obvious: 

(3) f (/ • 5u + m' ■ S(p — <5Il) dO + (p (F ■ 5u + M ■ 50) dl = 0. 

Here / and m' — force and moment loads per unit area (in dynamics they contain inertia 
terms), F and M — loads on the boundary contour dO, n — the density of energy of 
deformation. 
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For getting the stress tensors, deriving the equilibrium equations and Cauchy type formu¬ 
las, let’s use the technique, suggested in [12]. On virtual dislocations without deformation 


(4) 5n = 0, 5e = 0, 5^ = 0. 

After taking into account these equations in (3) we get variational problem with constraints 
(underlined). Introducing corresponding Lagrange multipliers the following equation arises 


( 5 ) 


I f ■ 5u + m' ■ 5(j) — t ■ 5e — n ■ 5x — A ■ (50 + V5n • n) 


dO 


+ <b (F-Su + M° • 50) dl = 0. 


Tensors r and /i are symmetric and belong to the tangent plane — such are constraints (4). 

Vector multiplier A (= Ay) is introduced according to (1) and allows to vary independently 
u and 0. 

Using the identities 

t • e = V • (l ■ u) ~ V ■ t ■ u 

fi ■ x = V ■ fi ■ 0 + n ■ b ■ u^j + V • n ■ 0-» V • • u 

A • Vu • n = V • (An • u) — V • (An) • u 
and theorem of divergence 


® v ■ udl = / (V ■ u T 2 Hn ■ u) dO 
J 80 Jo 

(v — ort, normal to dO in the tangent plane, 2 H = tr b — average curvature, u may be the 
tensor of any rank), let’s transform (5): 


( 6 ) 



(l + d ' h + A nj 


■ 5u 


m — ■ n j — A | • 50 


dO 


F_—v- + n ■ b + A nj | • 5u 


+ 1 M + v ■ d r' 


dl = o. 


Variations of 5u and 50 in the area and on the boundary are arbitrary, therefore the 
vectors in the curly brackets are zeros — that is the equilibrium equations and Cauchy type 
formulas (for F and M). This derivation relates not only to elastic behavior, term (—511) 
just denominates the virtual work of internal forces per unit area. 

But for elastic shell principle (3) allows to get the elasticity relations. Taking into account 
in (3) the consequences (6), we get: 


511 = r • 5e + /i • 5x 


=*► n 


n(f,^), 


T = 


an 

Tfe’ 




an 

dx' 










ON LAGRANGE THEORY OF SHELLS OF REVOLUTION 


4 


Elastic potential II in the linear theory is quadratic form, the number of coefficients in the 
case of general anisotropy is 21. In the simplest case of isotropy without cross couplings we 
get the relations with four constants 

2n = Cie 2 + C 2 e • e + Tfix 2 + D 2 x ■ x 
(7) e = tr e, x = tr x 

t = C\ea + C 2 e, fi = D \xa + D 2 h. 


In the limits of the above-said it is impossible to find the moduli C\-D 2 in elastic relations 
(7) — we should refer to asymptotic analysis of 3D problem by small thickness. But in the 
case of homogeneous isotropic material and independence of moduli from curvature we can 
take them like in the plate under plane stress and bending. 


Ci = 


2 nvh 
1-u 


h 2 

D x = —Ci 
12 


C 2 = 2 jih 


D 2 



2.3. System of equations and boundary conditions. Elasticity relations 

f ^' (= + E ' = + Qy) + / = o 

( 8 ) \ , x = " , 

(v ■ /fj + Q ~ w! — o 

(where Q = X — shear forces) and (7) resemble the known ones (about which, however, in 
the literature there is no common opinion), and coincide only with [2]. 

Boundary conditions, presented by the contour integral 


(9) 



V ■ 


Z + /i ■ h + Qn)) 


• 5u + (M° + v ■ fi) ■ 5(j) 


dl 


0, 


require reforming. Let’s introduce the tangent ort on contour dO — l = n x v. After partial 
integration we have 


| do i 

A ■ <pdl — ® [—A 


0| do = - {vd„u + Id (a) ■ n 
un ■ d v u + d^A ■ In) ■ u} dl. 


Transforming in this way the integral (9), we get the boundary conditions in the general form 


( 10 ) 


F — v ■ [l + ^ ■ b + Qn) + d l (A ■ In) 


■ 5u — A ■ un ■ d u 8u = 0 


where A = M + u ■ u. 

On the free edge the vector 5u and scalar n ■ d u 8u are unrestricted — their multipliers must 
be zeros. On the restrained edge u and n ■ d v u are defined — in all cases, including mixed 
ones, we have four conditions in components. 

It should be noted that during partial integration the contour was supposed to be smooth. 
Bit if there are corner points, then value d Z (A • In) contains ^-functions — like with point 
force in the load F. This well-known fact was being explained by many authors. 
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2.4. Consistency equations and statical-geometrical analogy. Components a a p and 
b a p must satisfy three differential equations of Gauss-Peterson-Kodacci [7] (following from 
symmetry r Q( g 7 by any pair of indexes). Their variation leads to the consistency equations 
for e and x. Let’s consider another derivation, connected with unambiguity of displacements 
and rotations. Let’s base on the theorem of circulation: 


( 11 ) 


dO 


(dr ■ v) dl — / n ■ V x vdO 
Jo 

v = Vn : n ■ V x u = 0 V • (n x (j) = 0. 

Let’s take the gradient of displacement in the form 


( 12 ) 


Vn = _e-axfl, 


Q = V x u ■ ( a + -nn ) = Q n n + n x 


By the reason of du = e ■ dr + D x dr, 0 is the vector of small rotation; but let’s note that 
generalized coordinates has only its “plain” part 0|| = U x = 0. Substitution of Vu from 
(12) in (11) leads to 

V • (n x e) — nV - fl + Vf2-n = 0 
(131 v- =/ - — — - 

Vf2 • n = —(V ■ (n x e))||. 

Next let’s consider “moment” tensor of deformation 


x = (VVn • 77 .) || = b ■ e — Vfi x n. 

Together with (13) it yields to 

VO = (x ~ b ■ e) x n — (V ■ (n x e))||n. 

The condition (11) for VO can be transformed to 

(14) V-(^*-6*-e* + An) = 0, (V • e*)|| - A = 0, 

where (...)* means the rotation of the tensor: x* = — n x x x n etc. 

Equations (14) express the consistency of deformation. A — is just a denotation, therefore 
we have three equations in components. Comparing (14) with (8), we see total analogy be¬ 
tween equilibrium and consistency equations. Statical-geometrical analogy is known already, 
but not in this form. 


3. Shells of revolution 

3.1. General equations. Let’s define the coordinate system on the shell: {0, s}, </> is the 
angular coordinate of meridional plane (longitude), s — arc coordinate along the meridian. 
Let’s express radius-vector in the form 

r = x(s)i + R(s)p((p). 

It should be noted that d^p = e^, = — p. Reference vectors in this coordinate system 

are: 

Zcj> = R(s)e,/, 

r s = x'(s)j + R'(s)p. 

Let’s pay attention to the following point: x and R are cartesian coordinates of meridian 
in the meridional plane, therefore ds 2 = dx 2 + dR 2 . Dividing both parts of this expression 
by ds 2 , we get 


x'(s) 2 + R'(s) 2 = 1. 
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The values x'(s) and R'(s) satisfy the trigonometric identity, that’s why the following deno¬ 
tation can be introduced: R'(s ) = sin ip(s), x'(s ) = cosip(s). The geometrical sense of ip(s) 
is the angle between tangent to meridian r s and symmetry axis i. Now reference vectors and 
normal may be written as 


= Re^ 

r s = i cos ip T p sin ip = t 
n = x t = p cos ip — i sin ip. 


Derivation formulas are 


d = — p = —t sin ip — n cos ip 
d^t = cos ipe^ 
d^n = sin ipe 4 , 


d s §<t> = 0 
d s t = ip'n 
d s n = —ip't. 


Introducing denotation R 1 = k, let’s build reciprocal basis and Hamilton operator 


r 4, = ke 0 
r s = t 

V = ke ^ + td s . 

Finally let’s calculate metric tensors of the shell of revolution 

a = Vr = / — nn = + tt 

b = -Vn = — k sin ipe^e^ + ip'tt. 

The equilibrium equations can be derived from the general ones easily. We’ll omit all 
further calculations. Here are only formulas for the gradient of vector and divergence of 
tensor. 

Vu = ( k (d^u^ + u s sin^ + u n cos^)) + tt ( d s u s — ip'u n ) 

+ §4 (k (d^u s - u$ sin ip )) + te# ( d s u^) 

+ e^n {k (d^Un - cos ip)) + tn ( d s u n + u s ip') 

V • p = [k {d'pH'f, + pfa sin ip + p s<t> sin ip) + d s p s<i> ) 

+ t(k (-/^ sin ip + fi s sin ip + d^pfa) + d s p s ) 

+ n {—kiitf, cos ip + n s ip'). 

On the base of these expressions and derivation formulas we can build the following se¬ 
quence of precise equations: 
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kinematic relations 


(15a) 


(15b) 


e<f, = k + u s sin ip + u n cos ip) 

e s = d s u s - u n ip' 

= e s <j> = ^ (k (d^Us - sin ip) + d s u^) 

X0 = k 2 (d‘ 2 )r pU n — 2cos ip — u s sin ip cos ip — u n cos 2 ip) 
+ k ( d s u n sin ip + u s ip' sin ip) 
k s = d 2 s u n + 2 d s u s ip' + u s ip" - u n ip" 

Ufa = k 2 (p—d'pUn sin ip + sin ip cos ip) 

+ k (i d 2 S(j) u n + d 4 ;u s 'ip' - d s u^ cos ip) 
x scj} = k ( d\ s u n - dgU^cosip + d^Ugip’) + k! {d <t> u n - u$ cosip) 


elasticity relations 


(15c) 

(15d) 


7"0 — Ci (e^ + e s ) + t s — Ci (e^ + e s ) + C 2 C S — r S 0 — 

= D\ (X(j> + x s ) + D 2 x r f ) Us — D\ (><0 + x s ) + D 2 x s 

/^0S X05 /^S0 k) 2 ^ S( j ) 


shear forces 


Qcj> = -k (< 9 0 /Z 0 + /J. 0.5 sin ip + /x S 0 sin - < 9 s /i S 0 + m s 
Qs = ~k (-/i 0 sin ip + /i s sin ip + - d s fi s - 


equilibrium equations 


(15f) 


" k 2 cos ip — fi S( f, cos ip sin ip) 

+ k (d^ + 2 t 0 S sin ip + /. i^ip ' sin ip - d s ii S(t) cos ip 

+ Us'pip' sin ip + cos + <9. S 70 S - cos ip + / 0 = 0 

/c 2 ^ cos ip sin ip + k (-70 sin ip + t s sin -0 + + fi s ip' sin ip + d^/i^ip') 

+ c^Ts + d s (i s ip' + n s ip" - Q s ip' + f s = 0 
k 2 ii ( j ) cos 2 ip + k (—70 cos ip + sin 

+ Tgip 1 + ii s ip' 2 + d s Q s + f n = 0. 


Let’s note that tensors x and /i are symmetric like r and e. It becomes obvious that X 0 S = 
x S 0 if we return from notation -0(s), k(s) to R(s), x(s). The equations for displacements 
can be derived from equations (15). Stress equations contain (15f) and consistency equations 
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(14), which in the case of shell of revolution are 


(16a) 


— /c 2 £ 0 S cos 0 sin 0 

+ k {d rj) x s - x s4> sin 0 - x^ sin 0 - <9 0 e s 0' + d s e^ s cos 0 + A^, cos 0 ) 

+ k'e^ cos 0 - <9, >Q, S = 0 
/c 2 £0 cos 0 sin 0 

+ k (-x s sin 0 + x^, sin 0 - d^x^ + e s 0 ' sin 0 + 

+ cos 0 — e.^0' sin 0 j + k's^ cos 0 + <9 S X0 — A s 0' = 0 
k (—x s cos 0 — £50^ cos 0 — £ < ^0 / cos 0 + c^A^ + A s sin 0) + X00' + <9 s A5 = 0, 


where 

(16b) 


A0 = /c (9 0 £5 - 2^5 sin 0) - a s £ 0s 

A s = k (-£5 sin 0 + £0 sin 0 - + < 9 S £ 0 , 


and strains are expressed through stresses by formulas 


(16c) 

(16d) 


x^ — Ai (p,0 + /o) + A2H<f> x s — Ai (/i^ + /i s ) + A2/i s 
X’^5 '^-2 f^(ps X00 

£0 = 5i (T0 + T s ) + -B 2 E0 £ s = B 1 (T0 + T s ) + B 2 T S £0 S = B 2 T^, S . 


Here 1?1 c 2 (2Ci+c 2 )- c 2 ’ d 2 ( 2 dI+d 2 )’ ^ 2 d 2 

Using formulas (15), given in this section, we’ll get the equations of cilyndrical shell, analyse 
them and compare by a number of properties with the equations of cilyndrical shells, given 
by other researchers ([3], [4], [14], [7], [ 6 ]) and with 3-dimensional model of the tube. 


3.2. Cylindrical shell. 


3.2.1. Equations. It’s not difficult to derive the equations of cylindrical shell from the equa¬ 
tions of shells of revolution (15). That’s enough to define the functions x(s ) and R(s) appro¬ 
priately. 


x(s ) — s = z, R(s ) = R. 

Function 0(s) can be calculated by formula 0(s) = arctan-^y = arctan^j. In this case 
it’s obvious, that 0 = 0 . 

Cylindrical shell is famous for its equations for displacements are differential equations with 
constant coefficients. Let’s change the variable for expression simplification s = (R , d s = -^d(. 
System (15) appears as 


(17a) 


££/(*C) + B(6C) = 0. 
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Here 


Aty(0, C)\ (u i(0,C)\ 

5(0,0 = u a ((j>,Q = w 2 (0, C) 
\u n {(j),c)j \u 3 i<f>,()j 

M,0\ (b i(0,C)\ 

(17b) 5(0,0= 0(0,0 = 0(0,0 

\O(0,O/ \O(0,O/ 



where 

Ln = + 2 y^ (d^ + 2<9| — L 12 = I /21 = — (1 + ^)d^d<j> 

^ ^ 5 22 = - (2<9^ + <9| — z/<9|) L 23 = L 32 = z/<9^ 

T 33 = 1 + 7 (AA + 1 — 2 <9| — 2z/<9|) L 31 = L 13 = <9^, + 27 (<9^ — <9^0, — <9|) 

O(0> C) = ^ (5/00,0 + "H(0,0) 

(17d) 0(0,0 = 5/00,0 

0(0, 0 = (5/n(0, 0 - 5 c m^(0, C) + <9 0 m s (0, 0) • 


Here z/, /1 — correspondingly Poisson’s ratio and shear module, A = (<9| + <9|) , 7 = 


3.2.2. .Roots of characteristic equation. Let’s consider sinusoidal loads and solutions 1 f/(0, 0 = 
U m ( 0e m *0 After substituting such solution in the system (17) and multiplying by e - " 1 * 9 / we 
get the system of differential equations with constant coefficients for displacement intensities 
f/ m (C). Trying complementary solution of the corresponding homogeneous system in the 
form U m (f) = U mk e k ^, is equivalent to trying the solution of input homogeneous system (17) 
in the form 

(18) [7(0, C) = U mk e kC+mi 0 

After canceling by e fc< ’ +m ®'0 we have the system of homogeneous algebraic equations for 
(U™ k , U™ k , . This system has solution only if its determinant equals to zero. 

During trying the solution in the form (18) differential operators turn into factors 

c 0 = im, c 0 = k. 

Thus resolvability requirement for such system can be written as 

V(im , k ) = 0, where 5(<9^, <9^) = det C. 


1 All the following is true also in the case of Fourier expansion of the solution U (c t>, C) = f2m =0 U m {C,)e m ' L< t 

with arbitrary loads. 
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Trying the solutions for U m (Q separately for each harmonic by 0, let’s represent this 
equation as equation for k. Coefficient m will be the parameter. 


(19) 


D = k 8 -7 — 2+ k 6 ( m 2 (2y — ( 1 / — 5 ) 7 2 ) + 1/7 + 4z/y 2 ) 

+ k 4 ^m 4 (—37 + 2{u — 2 ) 7 2 ) + m 2 (37 + 2(2 — v + ^ 2 )7 2 ) — 2^ 2 

4z / 2 — 5 1 \ 

+- 2 - 7 + 2 ^ ^ + ^ ) 

+ k 2 m 6 (27 — (y — 1)7 2 ) + m 4 (2(z/ — 1)7 2 — (4 + z/)y) 

+ m 2 (—(z/ — 1)7 2 + (2 + z^y) j — m 8 -7 + m 6 7 — m 4 -7 = 0. 


Let’s examine two cases: axisymmetric solution (m = 0) and antisymmetric 2 solution 
(m = 1). The solutions of the characteristic equation (19) for k in these cases are: 


m = 0 


^ 1 , 2 , 3,4 “ 0 ) ^ 5 , 6 , 7,8 “ 1 + 1 a / z ^ 2 T 


m = 1 


M, 2 , 3 , 4 — 0; ^ 5 ,6, 7 , 8 — 1 ^ ^2 + 5y + v + 3z/y ± 


I + 47 

A = — 7 ^ 1 — 5 — 47 + 97 2 + 4z^ + 22z/y + 3 OZ/ 7 2 
+ z / 2 7 - 1 + 9z/ 2 + 26z/ 2 7 + 25z/ 2 7 2 . 

Let’s look at the asymptotic form of these roots by 7 —» 0 (the shell is thin). 


k ^ -4 ± f 7 4 / 4 -y/ ±\/ z/ 2 — 1 + 7 1 / 4 


z/ 


/c 7 - 1 —> ± [ 7 1/4 v ±Vz / 2 — 1 + 7 174 


2 v / ±v / ^ 77 t 

2 + z/ 

2 v / ±v / ^ 77 t 


+ ... 


3.2.3. Solution. Let’s consider the case of sinusoidal load by both coordinates 

b m (<f>, z) = B m e ik ** +ik * z (m = 1,2,3 = 0, z, ra). 

Let’s consider also endlessly long cylinder, so that the border conditions on the ends may be 
neglected. The border conditions in the second coordinate 0 are periodical conditions and 
are satisfied automatically if we try the solution of the system (17) in the form 

u<j>((f>, z) = jj <l) e ik ’ l, ^ +ikzZ 

u z {<t>,z) = U z e ik ^ +ik * z 

u n (0, z) = ^k e ik ^ +ikzZ 
i 

(normal displacements have phase shift). 

2 Names “axisymmetric” and “antisymmetric” are quite conventional. Under additional conditions m = 0 

would mean either real axial symmetry, or torsion, and m = 1 should cause the bending of the cylinder like 

a beam. 
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In this case from (17) we get the system of linear algebraic equations 

r - £^2 (3 U^kl - 3U^R 4 k 2 z v + 6 U^R 2 k\ + h 2 U^R 2 k 2 z 

- h 2 U ( j ) R 2 klu + 2h 2 U ( f ) k0 + 3R 3 U z k c pk z + 3 R 3 U z k^k z u — 6U n k^R 2 — h 2 f/ n /c^, 

- h 2 U n k\ - h 2 U n k^R 2 k 2 z - 6B^R 2 ) = 0 

< - ^ (RU'pk'pkz - RU'frk'pkzis - 2U z R 2 k 2 z - U z k 2 , + + 2vRU n k z + 2B z ) = 0 

_ T2R^ UU ^ R2 ^ 2 ^^k(j> + 2 hkUffrkj, + 2h 2 U^k^R 2 k z 

+ Y2vR?U z k z - 12U n R 2 - h 2 U n - 2h 2 U n k\ - 2h 2 U n uR 2 k 2 z - h 2 U n R 4 k A z 

- h 2 U n k J - 2h 2 U n R 2 k 2 k 2 z - 12B n R 2 ) = 0 


In the case of normal load (B^ = 0, B z = 0, B n = — pR 0 \l — ) the solution for the 
displacement intensities U m: m = (j), z,n is 
(20) 

Utj, = ——j—^-(—HR 2 k^ + h 2 k^is — 12 i? 4 /c 2 — 2 h 2 R 2 k 2 z + QR 4 k 2 z u + QR 4 k 2 v 2 — 3 h 2 k 2 R 2 k 2 


+ Qk^R 2 u + h 2 k 2 ^v - 2h 2 R 4 k z + h 2 k 2 ^R 2 k 2 z u 


h 2 kl 


h 2 K) 


U z = ~ 6P ^ (~3h 2 k\v - QR 4 k 2 z v + QR A k 2 z v 2 + 6- Qkffiv - 2h 2 R 2 k\v + 2R 2 k 2 z v 2 h 2 
+ h 2 k 2 (j) + h 2 k^ + h 2 k^v + /r/c 2 .R 2 /c 2 + h 2 k 2 R 2 k 2 z v) 

U n = - ( ^-j(-2h 2 kl + 2h 2 klu + 6 k 4 z R 6 u - 6 k A z R e + 2 h 2 vR 4 kt - bh 2 klR 2 k 2 z + Uk^R 4 ^ 

- 12 k\R 4 k\ + 6 R?k\v - Gk^R 2 - h 2 v 2 R 2 k 2 z kl - 2h 2 R 4 k 4 z + 2h 2 k\R 2 k 2 z v) 

A = Qh 2 kl - 72kftRfv 2 - 24 k\R 4 h 2 v 2 + k 2 z h 4 k\ - 2k\h 4 k\ + k 2 z h% + 2Ak 2 z R 2 k\h 2 

— kl^k^u + 2 k z R 2 h 4 k^u — 2 k z R 2 h 4 k^v — Ak z R 2 h 4 k ^ + Ak 4 R 2 h 4 k 4 + 2Ak 6 z R & k^,h 2 
+ 5 k e z R 4 h 4 kl + 4 k\R 4 h 4 v - klR 4 h 4 k\v + 72 k\Rf + 6 k 8 z R 8 h 2 + 2k 8 z R 6 h 4 + 6 k 8 h 2 
+ 2 k 4 R 2 h 4 + 12h 2 k^R 2 k 2 v + 36k z R 4 k^h 2 — 2 k z R 2 h 4 k^iy 2 + 2k 2 h 4 k^ ) v — k 2 h 4 k^v 
+ 12 k\R?h 2 v + 30 h 2 R 4 k\ + 24h 2 ^i? 2 fc 2 - 12h 2 ^ 


- Y2h 2 R 2 k\k 2 z v - 36 h 2 R 4 k\k\ - 48 h 2 k\R 2 k 2 z 


3.3. Comparison with earlier theories. Equations of the shells of revolution, derived by 
us, differ from the results of other authors. Orders of equations in different variants are the 
same, and these equations are too complicated to draw conclusions about their qualitative 
resemblances and distinctions. That’s why let’s consider only cylindrical shell and some of 
its characteristics: those, that have been considered in the section 3.2. 


3.3.1. Operator coefficients. Operator coefficients of system of equations for cylindrical shell 
derived by different authors are given by C = £° + 7 C 1 . Major terms £° are the same: 
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L°n = \ {dl + 2 dl 

-^22 — 2 
l° 33 = i 

Smaller terms C 1 are: 


ud}) L ° v2 = L° 21 = ^( 1 + v)d c d </> 

v^l) l ° 2 3 = l ° 3 2 = 

rO _ r0 _ a 

-^31 — ^13 — °<t>- 


V. V. Eliseev [5] 

L \i = 2 (<9j + 2<9j - z/^) 

L 22 = 0 

L 33 = 1 + AA - 2dl - 2v&l L\ x = L\ 3 = 2 (1 - A) d 4> 

V. Z. Vlasov [14] 

L\ 1 = 0 L } 2 = L 21 = 0 

^22 = 0 L \3 = L 32 = ^ (<9|<9 C - 2 ^ - z/<9 2 <9 c ) 

-^33 = 1 + AA + 2c?0 I/ 31 = L 13 = — d^d^ 

A. L. Goldenveiser [7] & V. L. Biderman [3] 

L\i = 2d\ + d\ - 2 ^ I [ 2 = = 0 

L 22 = 0 L 23 = Lg 2 = 0 

-^33 = 2^A L 3l = L 13 = — ( 2d£ + d</, 

Yu. N. Novichkov [4] 


L \ 2 = Ll, = 0 


21 


-^23 — -^32 — 0 


L\ x = 2<9 2 + <9 2 - 2ud 2 c 
L 22 — 0 

-^33 = 


-^12 ~ — 0 
^23 = -^32 = 0 

^31 = L \ 3 = - 


(H + di 


2vdl) d 4> 


W. Fliigge [6] 

L\ x = |(i/ - 1)<9 2 L\ 2 = L\ l = 0 

£22 = \{v ~ l )$l l \ 3 = L k = ~ ( 2 <9c - dl + udl) <9 C 

-^33 = 1 + AA + 2<9^ L 31 = L 13 = - d'pdf- 


It’s hard to give preference to any of these formulas without regarding their grounds 
and peculiarities of their derivation. All variants of smaller terms have the same order 
of derivatives, similar level of complexity. All of them are symmetric. More significant 
differences will appear later. 
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3.3.2. Roots of characteristic equations. In the same way, as it was done in section 3.2 let’s 
examine the roots of characteristic equations in cases m — 0 and m = l 3 . 


V. V. Eliseev [5] 


m = 0 

^ 1 , 2 , 3,4 = 0 

4+7,8 = y/^T 1 + Wz/ 2 - 1 = 


± [ 7 1/4 v ±a/ z/ 2 — 1 + 7 1 / 4 


V 


2\J ±y/u'- 


2 -1 


m = 1 


zd 1 ) _ Q. 
"'1,2,3,4 — u > 


k 


( 1 ) 

5 , 6 , 7,8 


= ± 


l/lTn ( 2 + Y + i' + 3 ! / 7 ± Yi) 


± [ 7 ±V i/ 2 — 1 + 7 1 / 4 


2 + i/ 


2y/±VU- 


2 - 1 


+ 


A = —7 1 — 5 — 47 + 97 s + 4 z/ + 221/7 + 30 z/ 7 2 

+ i/ 2 7 _1 + 9i/ 2 + 26 z/ 2 7 + 25z/ 2 7 2 

V. Z. Vlasov [14] 

m — 0 

l(°) — n 

^1,2,3,4 — u 


(°) . - 4-'/ •/+ 


^5,6,7,8 — =*=7 


v/7 ± v 7-1 + z' 2 + 7 2 


z/ 


-1 + 7 

± ( 7 _ 1 ' / 4 \/±a/z / 2 — 1 — 7 1 / 4 


2 v / ±v / z^ : t 


+ 


m = 1 


4,2,3 — 0 


3 These characteristics are not necessarily derived in the referred books. We have derived these character¬ 

istics ourselves, using the equations, given in the books. 
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A. L. Goldenveiser [7] & V. L. Biderman [3] 


m = 0 

lS°) _ Q 
"T,2,3,4 — u 


k 


(o) _ 

5,6,7,8 

m = 1 


= ±7 ±Vv 2 — l 


uW 

,2,3,4 

U (!) 

'A,6,7,8 


= 0 

= ±7 -1 / 4 


\J rru ( 2v/ ^ ± v 7 - 1 + " 2 + 


± ( 7 V 4 V ±\/1/ 2 — 1 + 7 1 / 4 


Yu. N. Novichkov [4] 


m = 0 

u(°) _ Q 

^1,2,3,4 — u 

l(°) _ +^-i/4 

ft 5,6,7,8 — 




1 


m = 1 
^ 1 , 2,3 = 0 
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W. Fliigge [ 6 ] 


m — 0 


+°) _ n 

^ 1 , 2 ,3,4 ~ u 


k , 


(°) _ 
5 , 6 , 7,8 - 


y/l ± a/-1 + Z/ 2 + 7 2 


Z/ 


-1 + 7 

± ( 7~ 1 ^ 4 \/±V z / 2 — 1 — 7 1 / 4 


2V±v / z^ r t 


+ ... 


rn — 1 


l( 1 ) — n 

^1,2,3,4 — u 


k {1) - ±- 

"' 5 , 6 , 7,8 ~ 31 g 


7 (-1 - 27 + 37 2 ) 


± [ 7 1/,4 V/ ±V v 2 — 1 + 7 1 / 4 


A 7 ± 

2-i/ 


2v / ±7^ 


+ ... 


A = 8 + II7 + 97 s — 4 z/ — 15 z /7 — 9 z/y 2 
B = (A7) 2 — 47 (—2 — 47 + 67 s ) 

x (—2 — 87 — 67 s + 2 z/ 2 + 67Z/ 2 + 6 z/ 2 7 2 + 2 z/ 2 7 3 ) 


Here we can notice that when m = 0 all variants provide four zero roots for k. It means 
that corresponding complementary solution of homogeneous problem would be given by U — 
Co + CiC + C 2 C 2 + C 3 C 3 - Such polynomial solution resembles the solution of the problem of 
the Bernulli-Euler bending beam. Really, in axisymmetric case equations of cylindrical shell 
(17) are similar to the beam equations; each lengthwise fiber behaves like a beam. 

In the case rn — 1 there are significant differences in the solutions for k. Some variants 
provide beam-like solution, and some have only three zero roots. Beam solution in this case 
describes the deformation of the whole cylinder like a beam. Seems that such solution should 
exist. In some books tubes are frequently treated as rods [13]. Absence of such solution may 
be interpreted as substantial weakness of the system. 

Even more visible differences will be seen during comparing the solutions of the systems 
(17) by different authors with each other and with the solution of elasticity problem for 
3-dimensional tube. 


3.3.3. Solutions. Just like in the paragraph 3.2.3 we have built the solutions for U n under 
normal sinusoidal load in each of considered models. We have also examined the membrane 
(momentless) model that can be easily derived if we put 7 = 0 and built the solution for the 
3-dimensional cylindrical tube, loaded by the equivalent load on 

• internal surface: p in = pR/R in = pR/ (R— |) 

• outer surface: p out = pR/R ou t = pR/ (R + |) 

• both surfaces: p in = p out = pR/ (R in + R out ) = p/2, 

where pi n and p ou t are loads on correspondingly internal and external surfaces of 3D-cylinder, 
considered to be analogous to the only surface load in the shell. 

The solution of such 3-dimensional problem has the form u(cj),z,r) = [/(r)e* fc<#,9i+ * A;zZ for 
each component of the displacement vector. To compare the intensity of radial displacement 
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U r (r) with the normal displacement from shell models U n we take the averaged displacement 
TJ = - f R+h/2 u ( r ) dr 

Ur h JR - h /2 u r \' 1 Ul ' 

We have examined several cases of variability coefficients. They are shown on the figures 
1, 2, 3. 
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Figure 1 . Eliseev’s solution gives good estimate for displacements 
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FIGURE 2. Inner load differs from the both other 

On these figures the dependence of the product of normal displacements and the thickness 
U n (h)h from thickness h is shown. We choose U n (h)h instead of U n (h ) because the function 
U n (h ) is hyperbolic, has singularity in h = 0 and spans across too big range in the ordinate 
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1.2 
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Figure 3. All but two shell models give very close solutions 


axis, making the details of function behavior too small on the plot. Multiplying by h removes 
hyperbolic aspect of the function and makes the plot more comprehensible. 

So, first of all we may check that in all cases small thicknesses make the results of the shell 
and membrane models the same. That is natural, because when 7 —> 0 all shell models lose 
their personal £} terms and turn in the membrane model indeed. 

But 3D-models give a bit higher values of radial displacement in some range of low h - 
rigidity occurs to be lower than in 2D-models. This may probably be explained by the less 
number of modeling suppositions, conditions, restraints on the nature of the deformation 
process in the classic 3D elasticity. 

Next, we may see that membrane model always gives higher displacements than most of 
the shell models and 3D-models in higher range of h. That is caused by the higher rigidity 
of 3D and shell models, brought by bending (moment) terms. 

In all cases we may also notice that 3D-cylinder with inner load differs from the both other 
3D-models, especially in the higher range of h. We can conclude that the model with inner 
load seems to be somewhat special, because it leads to the singular load in limit case h = 2 R. 

In the axissymmetric case (k^ = 0) we may see that Eliseev’s solution stands apart from 
all other shell models, but it gives good estimate for displacements also in the very high 
h-range, reflecting the possibility of so special inner load case. 

In the antisymmetric case (k# = 1) all but two shell models give very close solutions. 
“Outsiders” are Vlasov’s and Fliigge’s models — they are discontinuous. Let’s examine this 
discontinuity. 

The expressions for the U n in both these models are rational functions as well as (20). And 
due to the polynomial in h in the denominator they may have discontinuities. In all models 
there is multiplier h in the denominator of U n , so this function is discontinuous in h — 0 . 
But there is also another point of discontinuity. Using MAPLE it is quite easy to find them. 
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In the case k$ = 1, R = 1, v = 0.25 these points are given by the equations: 

Vlasov — h d : (14*4 + 8/c® + 19*4 + 3 )h d 

+ (144 - 96 k 6 z - 96 k\ - 336 kj)h 2 d - 1080 k\ = 0 
Fliigge — h d : (-18*4 + 73*4 + 48*:®)/i® 

+ (-1212*4 - 1608*4 - 768 k\ + 288 )h\ 

+ (-6912 + 8064*4 + 2304*:® - A176k 2 z )h 2 d + 25920*4 = 0 

Though Vlasov’s model gives biquadratic equation for h d with taken parameters, in general 
case there’s no way to solve these equations analytically, therefore we may find the real roots 
numerically. On the figure 4 there is the dependence of h d , where the discontinuity takes 
place, from k z . It should be noted also that the second elastic module p has no influence on 
this relation. 



Figure 4 . Discontinuity thickness dependent on k z 

We can see that discontinuity takes place in the possible range of h e (0; 2 R) and corre¬ 
sponds with figures 2, 3. 


4. Conclusion 

In this study we have faced with high correlation of shell displacements with 3D ones 
also in thick shells, where any simplifying preconditioning numerical algorithms (see e.g. 
[10]) would not be useful, because of the acceptable condition number, and computational 
resources would be wasted without big improvement in precision. 

On the other hand such approach to the shell theory doesn’t give an idea about the 
connection of the stress tensors r and p with 3D stress, therefore we cannot use these solutions 
within the well-developed strength conditions. This is the matter of the future studies as 
well as the derivation of other special cases of shell theory, e.g. shallow and narrow shells. 

It is obvious that direct tensor algebra makes the analytical calculations much easier, 
bringing besides research advantages the pedagogical ones also. 
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The surprisingly inaccurate results, given by some of the considered shell models, make us 

understand that the shell theory is not well formed and closed even in its basis yet. Perhaps, 

even more unsunned surprises will be discovered during the foregoing future studies. 
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